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Abstract 

When taking the real, inhomogeneous and anisotropic matter dis- 
tribution in the semi-local universe into account, there may be no need 
to postulate an accelerating expansion of the universe despite recent 
type la supernova data. Local curvatures must be integrated (over 
all space) to obtain the global curvature of the universe, which seems 
to be very close to zero from cosmic microwave background data. As 
gravitational structure formation creates bound regions of positive 
curvature, the regions in between become negatively curved in order 
to comply with a vanishing global curvature. The actual dynamics of 
the universe is altered due to the self-induced inhomogeneities, again 
more prominently so as structure formation progresses. Furthermore, 
this negative curvature will increase as a function of time as struc- 
ture formation proceeds, which mimics the effect of "dark energy" 
with negative pressure. Hence, the "acceleration" may be merely a 
mirage. We make a qualitative and semi-quantitative analysis, using 
newtonian gravity corrected for special relativistic effects, which works 
surprisingly well, to corroborate and illustrate/ visualize these state- 
ments. This article may be seen as a plea to start taking seriously the 
observed inhomogeneous distribution and the nonlinearities of non- 
perturbative general relativity, and their impact on the dynamics and 
behavior of the cosmos. 
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Measurements since the late 1990s on type la supernovae (SN) PP,[2] sur- 
prisingly seemed to indicate that the universe accelerates its expansion at 
the present epoch, instead of the deceleration expected if gravity is univer- 
sally attractive. This finding has resulted in a "neo- standard" interpretation 
where the present universe is believed to be dominated by "dark energy" with 
negative pressure. However, when taking the observed inhomogeneous struc- 
ture of the universe into account, and considering the real geodesic paths of 
observed SN light, such a hypothesis might be superfluous. 

i) If we assume that luminous matter (and by necessity gas, dust and 
plasma for star formation) is a good "tracer" of regions of higher than aver- 
age density, almost all photons from distant objects that reach us on earth 
must have traversed regions with little or no matter with which it can in- 
teract electromagnetically. The light from distant objects "zigzag" through 
the maze defined by gravitationally bound objects (galaxies, galaxy clusters, 
superclusters) making the real (geodesic) path longer than the one calcu- 
lated from standard, perfectly isotropic and homogeneous, FRW-cosmology, 
the more so the longer structure formation has progressed (z <~ 1) [3]. ii) 
The actual dynamics of the universe is altered due to the self-induced inho- 
mogeneities, again more prominently so as structure formation progress. We 
will, deliberately from a pedagogical standpoint, employ a heuristic model 
using newtonian gravity, with special relativistic corrections, both for the 
ease of visualization/interpretation and the much simpler (linear) mathe- 
matics than full-blown (nonlinear) general relativity, in which not even the 
two-body case is analytically solvable. Even though newtonian gravity is 
not completely mathematically consistent in an infinite universe, due to its 
dependence on boundary conditions infinitely far away, we can side-step this 
by treating the Big Bang as an "explosion" in a pre-existing (newtonian) 
space. As can be seen in Fig. 2, this special relativistically corrected newto- 
nian model comes surprisingly close to the general relativistic FRW-model, 
especially for low z as expected. 

In general relativity, gravitationally bound systems have a positive space- 
time curvature. At the same time we know, from observations of the cos- 
mological microwave background radiation (CMBR)[1],[S], that the global 
geometry of the universe most probably is fiat. This means that the curva- 
ture between gravitationally bound systems (solar systems, galaxies, galaxy 
clusters, etc) must be negative. This conclusion applies to all globally flat 
universes with (semi-)localized gravitationally bound systems. 

For a truly exact description, we would need to know the energy-momentum 
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tensor (T^u) at each point between us and the distant SN, which is physically 
impossible. And even if we had such perfect information, it would still be 
mathematically inconceivable to solve the resulting Einstein equations, 



(1) 



to deduce the local curvature at each point in terms of the Riemann curvature 
tensor, due to the complexity of the equations - ten coupled, nonlinear PDEs. 
(The real lure of assuming perfect homogeneity and isotropy is that Einstein's 
eqs. simplifies to two coupled linear ODEs, Friedmann's eqs., in which the 
total dynamical behavior of such FRW-model universes is contained solely in 
a(t), the cosmic scale factor.) 

One hope would be that an approximate model could be used to deter- 
mine an "effective" curvature for bound systems and for the space in-between. 
It would define a semi-local mean curvature parameter, fc, for the regions 
bound/between. (This parameter is related to the scalar curvature R = R^, 
averaged over the region, {R).) Even in such a "coarse-grained" inhomoge- 
neous model, the cosmic scale factor a{t) will be scale-dependent, as noted 



In terms of = p/ pcru (where pcrit is the density required for flatness): 



Or, stated in terms of the mass-energy density: p global = Pcrit, Pbound > Pcru, 
Pbetween ^ Pcrit- Howcvcr, ouc should keep in mind that the very definition 
and usage of fl assumes homogeneity and isotropy. 

As matter preferentially clumps in well localized objects (stars, etc), the 
majority of the photons that reach us travel mainly in "under-dense" (neg- 
atively curved) f2 < 1 space. By observing light we are thus automatically 
biased to measure an "apparent" curvature which is less than the actual 
global curvature|§ Light from a SN source will, in a universe which exhibits 

^ "This leads one to view the Hubble parameter as a stochastic variable, subject in the 
hierarchical scheme to effects of local density fluctuations on all scales." 

^Neutrinos and gravitational waves should show less bias as they can travel unhindered 
through huge amounts of matter without interacting appreciably. This means that both 
negative and positive curvature effects contribute, which reduces the bias. 





(2) 



^bound ^ 1) 
'between ^ 1 



(3) 
(4) 
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gravitational clumping/structure formation, always be switched towards a 
seemingly more negatively curved universe. For a universe with zero global 
curvature, the SN light will thus approach the curve for an open universe, 
see Fig. 1. 

Another compelling property is that this negative curvature effect will 
automatically mimic a very small cosmological constant, beginning to "dom- 
inate" at an epoch when a significant amount of structure has evolved. Be- 
fore structure formation through gravitational condensation becomes effec- 
tive (1100 3> -2 > 4), all space will have roughly the same curvature (A; ~ 0). 
However, structure formation will produce bound systems with increasing 
^bound-i which means that ^between will be a decreasing function of time. 
Hence, the space in-between bound systems will asymptotically approach 
n = with time (under-density being diluted by expansion), simulating an 
accelerated expansion. 

After these general considerations, we now turn to the newtonian model 
to qualitatively and semi-quantitatively, make our case about ii), the dy- 
namics. First, it is possible to deduce analytical results from a maximally 
inhomogeneous and anisotropic distribution - a two-body problem: 

Nomenclature: 

(i, the distance from an object to the observer. 

the observed velocity of an object in the radial direction, 
a, the observed acceleration of an object in the radial direction. 

We first define the "homogeneous acceleration", denoted ah- This is the 
acceleration that would be experienced by an object if the distribution would 
have been homogeneous and isotropic on all scales. Simply put, if a test par- 
ticle is somewhere in a homogeneous sphere, the only net effect of gravitation 
is the mass within a smaller sphere with a radius equal to the distance from 
the test mass to the center. The acceleration for this test mass becomes 

= (5) 

where m is the mass of the small sphere. This mass can be expressed by the 
volume ratio times the total mass of the distibution, M, 

m = M^, (6) 
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Figure 1: Shown [7] are the original data points of the High-z Supernova team 
(filled squares) in, and the Supernova Cosmology Project (open squares) [2]. The 
dashed line is the theoretical prediction for a homogeneous and isotropic (FRW) 
universe which is flat and without cosmological constant (Om = 1) ^^A = 0). The 
solid line is the corresponding prediction for an empty (open) universe {^Im = 0, 
= 0). Also shown (short dashes) is the theoretical prediction for a flat universe 
with solely a dark energy component (17 Af = 0, Oa = !)• The dotted line is the 
solution currently favored for the SN la data by both experimental groups (Qm = 
0.3, = 0.7). For z > 1, where unfortunately also observational measurements 
become increasingly difficult, it starts to deviate towards the {Qm = 1, ^^A = 0) 
line. The neo-standard explanation for this is that f^A has become dominant only 
fairly recently {z < 1). Observed SN photons within a globally flat universe will 
always tend towards the line for the open universe, due to inhomogeneous structure 
formation (assuming I^a = 0). 
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where it! is the radius of the total volume. The acceleration can be expressed 
as 

GMd 

= (7) 

If we know the total mass, then one way of deciding R is to check the mean- 
distance of the observed supernovas. If we assume that we can see all super- 
novas in the observable universe, we can then just take the mean-value of d 
and multiply with four-thirds to get the radius, as the geometric center lies 
at three-fourths of the radius in a cone-fragment of a homogeneous sphere. 
The equation for the acceleration becomes 

27GMd 

where (d) is the observed mean-value of the distance, d. 

Another way is if we know the mean-density, p. Then the acceleration 
becomes 

ah = — • (9) 

Let us now obtain a quantitative measure of inhomogeneity. If we start 
with the acceleration, we can simply add and subtract the "homogeneous 
acceleration", ah, 

a^ah + {a- Qh) , (10) 



2 

extract a factor ^ from the bracket 



fad ahd\ , , 



Introducing the Hubble parameter, if = ^, we get 

a^a,^m^[^-f)^ (12) 

The terms inside the bracket are dimensionless. As can be seen, this inho- 
mogencous result completely without dark energy yields the same behavior 
as the standard (homogeneous) model with a cosmological constant, 

a^ah + dH^ftA. (13) 
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where instead of the bracket one has the Q/i^ term ^U\. From now on the 
correction term for the inhomogeneity will be denoted by 

Let us consider an inhomogeneous Universe consisting of two large bodies 
and an observer, whose mass is negligible, situated between them. The two 
bodies have equal mass and therefore half the mass of the Universe. 

As this is a two-body problem the acceleration in the radial direction is 

(15) 



The gravitational potential, V, is defined according to 

F = -S72d V, (16) 
where F is the force acting on the particles in the radial direction, 



16^2 • 

Solving for V gives 



(17) 



— 

If the particles have escape velocity (corresponding to flat space-time), the 
virial theorem states that 

2K + V = 0, (19) 

where K is the kinetic energy of one of the particles, K = Solving for 

f 2 gives 

9 GM , , 

The total mass of the Universe is M so, 

27GMd 

In this case (rf) = d because both particles are at the same distance from the 
observer. So the "homogeneous acceleration" becomes 

ay, = (22) 
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1 27\ 

Q={ — + — = 1.1875. (23) 



which gives the value, 



2 16/ 

Thus the real acceleration will be 

a = ah + dHh.l875. (24) 

The motivation for this somewhat trivial demonstration is to show that if 
we observe an inhomogeneous Universe and still "pretend" it is homogeneous, 
it will yield a correction term for the inhomogeneity that behaves in the 
same way as a cosmological constant. This points to a possibility to remove 
the need for a mysterious dark energy which today has no fundamental, 
microscopic explanation or justification whatsoever. 

An N-body simulation with newtonian gravitation that takes into account 
the effects of special relativity has also been done. The particles have a small 
gas-like extension so if they get really close to each other the gravitational 
potential will not be infinite. Instead of having the ordinary potential going 
like — i (where d now is the distance between two particles), it will go like 
if the particles are closer than their radius, r, given at the start of the 
simulation. If the particles are separated by a distance larger than r there 
will be a potential like 

V = -G^, (25) 
where m is the mass of each particle, while for < r it is 

2r r 



V^-G^ + G-^. (26) 



The particles will therefore not collide but simply pass through each other 
without having an infinite acceleration, just as expected for the "test-particles" 
of cosmology; galaxies. They are distributed randomly in a sphere so that 
the large-scale density is almost uniform. The particles begin with escape 
velocity as initial condition, as this is equivalent to a flat Universe in the Ein- 
stein setting. Then the acceleration is calculated and the next distance and 
velocity for all particles are iterated using a finite difference method. Plotting 
the distance from the center of the sphere against red-shift, and comparing 
to a few general rclativistic "standard" cosmological models, show some in- 
teresting results. Fig. 2. All models give similar behavior up to z ~ 1 and 
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Figure 2: The prediction of three different homogeneous models are shown to- 
gether with our simulated inhomogeneous "particles". The solid line is an open 
FRW-model with = 0.2, = 0. The dotted hue is a flat FRW-model with 
Oftf = 1, = 0. The dashed line is the (flat) newtonian prediction with cor- 
rections for special relativity. An inhomogeneous matter distribution may thus 
be interpreted wrongly as to suggest that VIm is lower than it actually is, e.g. 

= 0.2 if interpreted through a homogeneous and isotropic "standard model" 
cosmology. The simulated particles correspond to Q = Qm = ^- The assump- 
tion of homogeneity and isotropy mislead the physical interpretation if the real 
distribution is inhomogeneous, as is the case for our Universe. 
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the particles in the flat but inhomogeneous model approach the homogeneous 
open model for higher z, in line with our previous analytical result. 

We end with some related comments: 

• SN la data probe regions with z < 1.7. Homogeneity and isotropy 
is valid only on scales significantly larger (orders of magnitude) than 
the cosmological "voids" and "filaments" [HIIS], i-^-, at distances 3> 120 
Mpc (corresponding roughly to z ~ 0.03). Neither can one a priori rule 
out clumping on even grander scales. Hence, the "cosmological prin- 
ciple" of homogeneity and isotropy, which the FRW-solution crucially 
depends upon, does not apply exactly. Instead full consideration of the 
inhomogeneities should be taken, at least up to the distance scale where 
homogeneity and isotropy may be considered a valid approximation. 

• The CMBR almost certainly probes the overall geometry/curvature of 
the universe [z ~ 1100), as little gravitational structure could form/grow 
before photon decoupling. The statistical weight of the low 2;-range 
where appreciable structure has formed is negligible compared to the 
higher z-range which thus dominates the integrated effect for the CMBR. 
For very high redshift the photons accordingly should behave "as ex- 
pected" in a fiat universe. Also, the CMBR is "everywhere" while SN 
photons travel from a pointlike source to us along a sharp geodesic 
"ray". This means that, due to the inhomogeneity and anisotropy 
at small to medium scales, constraints from SN and CMBR may not 
"carry over" trivially between one another. 

• In a FRW-universe with a cosmological constant it is just a strange 
and completely unexplained "cosmic coincidence" that Qm ^ now 
{Qm ^ earlier and Qm ^ ^^a later). However, in our scenario it is 
an automatic bonus, as an appreciable amount of structure must form 
before intelligent life can evolve to observe it. It is thus natural that 
we live in an epoch when the apparent "acceleration" (really due to 
inhomogeneity) becomes observable. 

In conclusion, we have noted that by regarding the real inhomogeneous 
matter distribution arising from time- dependent gravitational structure for- 
mation, it might be possible to avoid the conclusion that the expansion of 
the universe accelerates, as normally drawn from high-z SN la data. This 
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would alleviate the need to postulate that the present universe at large is 
dominated by an exotic "dark energy" with a mysterious negative pressure. 

This can also be seen as a plea to start taking seriously the nonlinear 
character of general relativity, one of the intrinsically most highly nonlinear 
theories in existence, where even very "small" perturbations could grow ex- 
ponentially. This, and other phenomena, are known to occur in other (e.g. 
chaotic) nonlinear systems, whereas overly simplified linearized models, or 
even the whole perturbative series with infinitely many terms, suppress or 
even exclude such behavior. Coupled with the vast distances and timescales 
relevant in cosmology this could make all the difference. 

For example, a FRW-model with newtonian perturbations as used in sim- 
ulations of large-scale structure formation may be overly simplified to capture 
the true dynamics. We believe that the nonlinear aspects of gravity have been 
gravely underestimated in cosmology, including very important but unsolved 
turbulent processes. 

Despite high-precision observations advancing cosmology towards becom- 
ing a normal "exact science" in the last couple of decades the present physical 
understanding of the universe as a whole is probably still very crude. We 
predict that theoretical cosmology fifty years from now will bear little resem- 
blance to the FRW-models almost universally used today. 
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